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HIGHLIGHTS 


►  Particle  size  and  volume  fraction  in  SOFC  electrodes  strongly  influence  performance. 

►  Binary  mixtures  of  spherical  electron  and  ion  conducting  particles  in  SOFC  analyzed. 

►  We  propose  a  new  estimate  for  the  three-phase  boundary  reaction  domain. 

►  Parameter  study  of  morphology  for  enhancement  of  electrode  performance. 

►  Unlike  sized  particles  in  cathode  increase  SOFC  performance  up  to  15%. 
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Optimization  of  the  microstructure  of  porous  electrodes  plays  an  important  role  in  the  enhancement  of 
the  performance  of  solid  oxide  fuel  cells.  For  this,  microstructural  models  based  on  percolation  theory 
have  proven  useful  for  the  estimation  of  the  effective  material  properties  of  the  electrode  material, 
assumed  to  consist  of  a  binary  mixture  of  spherical  electron  and  ion  conducting  particles.  In  this  work, 
we  propose  an  extension  of  prior  approaches  for  calculating  the  effective  size  of  the  three-phase 
boundary,  which  we  judge  to  be  physically  more  sound  and,  in  particular,  well  suited  for  character¬ 
izing  mixtures  of  particles  of  different  sizes.  This  approach  is  then  employed  in  a  one-dimensional  cell 
level  model  encompassing  the  entire  set  of  processes  of  gas  transport,  electronic  and  ionic  conduction  as 
well  as  the  electrochemical  reactions.  The  impact  of  the  electron  and  ion  conducting  particle  sizes,  their 
volume  fraction  and  their  size  ratio  on  the  performance  of  the  fuel  cell  are  investigated  in  a  parametric 
study.  Under  certain  conditions,  cathode  microstructures  having  electronic  conducting  particles  of  size 
different  from  that  of  the  ionic  conducting  particles  become  preferable  and  yield  a  higher  maximum 
power  density  when  compared  to  the  best  possible  configuration  of  monodisperse  particles. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cells  (SOFC)  directly  transform  chemical  energy 
from  gasified  fuels  into  electrical  and  heat  energy  in  a  highly  effi¬ 
cient  way.  The  microstructure  of  the  composite  anode  and  cathode 
plays  a  critical  role  in  the  high  performance  of  the  SOFC.  Consisting 
of  binary  mixtures  of  electron  and  ion  conducting  particles  that 
form  a  porous  composite,  the  electrodes  must  sustain  the  transport 
of  electrons,  ions  and  the  gaseous  reactants  and  products  from  and 
to  the  electrocatalytic  active  reaction  sites,  i.e.  a  small  region  close 
to  the  interface  between  percolated  electronic  and  ionic  particles 
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and  exposed  to  the  pore  space,  known  as  the  three-phase  boundary 
(TPB). 

A  SOFC’s  performance  is  limited  by  different  loss  mechanisms. 
Activation  losses  occur  for  the  hydrogen  oxidation  reaction  in  the 
anode  and  the  oxygen  reduction  reaction  in  the  cathode,  when 
electrons  are  transferred  between  the  respective  conduction  bands 
and  gas  molecular  orbitals  at  the  TPB  sites.  Ohmic  losses  arise 
predominantly  in  the  transport  of  oxygen  ions  through  the  elec¬ 
trolyte  from  the  cathode  to  the  anode  reaction  sites,  but  also  for 
electron  transport  between  the  current  collectors  and  the  TPB 
through  the  percolated  electron  conducting  particle  network. 
Finally,  the  diffusive  gaseous  transport  of  oxygen  through  the 
porous  cathode  as  well  as  hydrogen  and  water  vapor  through  the 
porous  anode  is  governed  by  concentration  losses. 
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These  loss  mechanisms  can  be  directly  correlated  to  the 
morphological  properties  of  the  microstructure.  Activation  losses 
are  strongly  influenced  by  the  disposable  reaction  sites,  i.e.  the 
amount  of  TPB  in  the  electrodes,  which  in  turn  is  linked  to  the 
respective  electronic  and  ionic  conducting  particle  sizes  and 
volume  fractions.  Ohmic  losses  are  governed  by  the  effective  elec¬ 
tronic  and  ionic  resistivites  of  the  porous  electrode  material  that 
depend  strongly  on  the  volume  fractions  and  percolation  proba¬ 
bilities  of  the  respective  particle  networks.  The  pore  radii  and 
porosity,  also  depending  on  the  particle  sizes  and  volume  fractions 
in  the  porous  electrodes,  directly  affect  the  gas  transport  from  the 
fuel  channels  to  the  TBP  and  consequently  the  concentration  losses. 
Hence,  by  optimization  of  the  morphological  parameters  of  the 
microstructure,  i.e.  the  diameters  of  the  electronic  and  ionic  con¬ 
ducting  particles,  their  size  ratios  and  their  volume  fractions,  the 
losses  can  be  controlled  and  the  performance  of  the  SOFC  can  be 
improved. 

For  designing  highly  efficient  fuel  cell  electrodes,  the  correlation 
between  the  parameters  of  their  microstructures  and  the  resulting 
disposable  active  TPB  for  the  electrochemical  redox  reactions  is  of 
special  interest.  Numerical  modeling  techniques  enable  us  to  gain 
insight  into  this  correlation:  One  possible  approach  is  to  numeri¬ 
cally  generate  electrode  microstructures  by  three-dimensional 
ordered  or  disordered  packing  of  spherical  particles  and  subse¬ 
quently  calculate  all  relevant  parameters  from  these  structures 
[1—4].  In  another  class  of  model,  which  we  employ  in  this  work,  the 
SOFC  microstructure  is  described  based  on  an  assumed  average 
particle  coordination  number  and  percolation  theory  [5-11  ].  While 
the  3D  particle  packing  approach  does  not  need  an  assumption  for 
a  constant  average  coordination  number,  it  requires  extensive 
computational  efforts  compared  to  the  approach  based  on  perco¬ 
lation  theory.  Thus,  the  latter  becomes  especially  favorable  when 
examining  different  microstructures  having  a  wide  range  of  particle 
sizes  and  volume  fractions. 

Two  different  formulations  of  the  effective  size  of  the  TPB  can  be 
found  in  literature.  They  differ  in  the  assumption  regarding  where 
the  electrochemical  reaction  takes  place.  Costamagna  et  al.  [5] 
introduced  the  TPB  as  an  area  per  unit  volume  (m2/m3)  param¬ 
eter,  being  proportional  to  the  contact  area  of  intersecting  elec¬ 
tronic  and  ionic  conducting  particles.  While  this  first  percolation 
micro-model  for  SOFCs  yields  realistic  estimates  for  the  TPB  size 
that  are  successfully  employed  in  a  number  of  theoretical  studies 
[6,8,10—13],  it  is  based  on  the  assumption  that  the  electrochemical 
reaction  happens  on  the  interparticle  contact  area  that  is  not  con¬ 
nected  to  the  pore  phase  and  therefore  inaccessible  to  the  gaseous 
reactants.  In  contrast,  the  TPB  has  also  been  formulated  by  Sunde 
[14]  and  others  [1,7,9,15,16]  as  a  length  per  unit  volume  (m/m3) 
parameter,  hence  suggesting  that  the  electrochemical  reaction 
takes  place  on  a  (one-dimensional)  circular  contact  line  at  the 
intersection  of  the  surfaces  of  electronic  and  ionic  particles.  This 
approach  has  been  shown  to  be  in  good  agreement  with  experi¬ 
ments,  i.e.  measurements  of  the  TPB  length  from  focussed  ion- 
beam  scanning  electron  microscopy  imaging  the  three- 
dimensional  microstructure  [17—20],  It  also  predicts  TPB  scaling 
consistent  with  experiments,  i.e.  an  increase  of  the  TPB  by 
approximately  a  factor  of  four  when  the  particle  diameters  are 
reduced  by  half  [21—23],  The  drawback  of  this  approach  is, 
however,  the  assumption  that  the  electrochemical  reaction  site  is 
hypothetically  confined  to  a  one-dimensional  line  at  the  interpar¬ 
ticle  junction.  Various  recent  studies  based  on  density  functional 
theory  [24]  and  models  coupling  elementary  chemistry,  electro¬ 
chemical  oxidation  reaction  kinetics  and  diffusion  mass  transport 
phenomena  [25-27]  indicate  that  the  electrochemical  reaction 
rather  takes  place  on  the  particle  surfaces  within  a  finite  distance 
from  the  contact  line.  In  order  to  accommodate  these  findings,  we 


propose  an  extension  of  the  percolation  micro-model  that  is 
physically  sounder  by  allowing  the  electrochemical  reaction  to  take 
place  on  the  particle  surfaces  within  a  finite  width  of  the  inter¬ 
section  between  electronic  and  ionic  particles. 

The  extended  percolation  micro-model  is  then  utilized  as 
a  foundation  for  calculating  effective  material  parameters  for  the 
porous  anode  and  cathode,  such  as  electronic  and  ionic  conduc¬ 
tivities,  active  three-phase  boundary  area  and  mean  hydraulic  pore 
radii.  Since  a  main  objective  of  this  paper  is  identifying  morpho¬ 
logical  microstructural  parameters  that  enhance  the  performance 
of  a  SOFC  system,  an  ideal  balance  among  the  ohmic,  activation  and 
concentration  losses  is  sought.  Hence,  a  one-dimensional  cell  level 
model  is  employed  to  systematically  investigate  the  impact  of  the 
electrode  microstructure  parameters,  such  as  the  electronic  and 
ionic  conducting  particle  diameters  and  their  volume  fractions,  on 
the  overall  performance  of  the  fuel  cell  system.  This  model  is 
parameterized  through  experimental  data  of  Jiang  and  Virkar  [28], 

Multiple  theoretical  studies  showing  the  effect  of  electrode 
microstructure  on  losses  in  solid  oxide  fuel  cells  can  be  found  in 
literature:  Chan  and  Xia  [6]  as  well  as  Cannarozzo  et  al.  [13] 
investigated  the  influence  of  anode  particle  diameter  on  activa¬ 
tion,  ohmic  and  concentration  losses,  assuming  equal  sized  elec¬ 
tronic  and  ionic  particles  and  the  effective  material  parameters  of 
Costamagna  et  al.  [5].  In  their  investigation  of  the  overpotentials  of 
cathode  interlayers,  Chen  et  al.  [17]  considered  the  effect  of 
different  particle  diameters,  size  ratios  and  volume  fractions.  In 
addition,  there  are  various  parameter  studies  for  microstructure 
optimization  of  both  conventional  [7,10,29]  and  functionally  graded 
[8,11  ]  electrodes.  The  main  focus  in  these  studies  is  on  finding  the 
ideal  thicknesses  of  anode,  dense  electrolyte  layer  and  cathode  and 
the  ideal  particle  diameter,  while  relatively  little  attention  is  paid  to 
the  effect  of  the  volume  fractions  and  the  particle  size  ratio  of  the 
electronic  and  ionic  conducting  phases  on  SOFC  performance 
[30,31].  With  only  a  few  exceptions,  losses  are  calculated  for  only 
a  single  constant  operating  current  density,  and  the  effects  of 
specific  microstructural  properties  on  losses  are  examined  in  an 
isolated  manner  for  either  the  anode,  electrolyte  or  cathode  rather 
than  addressing  the  microstructure’s  impact  on  the  performance  of 
the  fuel  cell  as  a  whole. 

Hence,  the  main  objective  of  this  paper  is  to  systematically 
investigate  the  influence  of  particle  diameter,  size  ratio  and  volume 
fraction  in  anode  and  cathode  on  effective  material  parameters 
such  as  TPB,  electronic  and  ionic  conductivities  and  hydraulic  pore 
radii  and  on  the  performance  of  the  fuel  cell  as  a  whole.  Employing 
our  new  approach  for  estimating  the  extent  of  the  TPB,  the  impact 
of  morphological  parameters  on  the  fuel  cell’s  performance  will  be 
demonstrated,  and  strategies  for  enhancing  the  performance  will 
be  presented. 

2.  Method  of  approach 

2.1.  Description  of  anodic  and  cathodic  redox  reactions 

In  this  study,  we  consider  an  anode-supported  SOFC  design  with 
single-layer  electrodes  as  shown  in  Fig.  1.  It  consists  of  a  nickel- 
yttria  stabilized  zirconia  (Ni-YSZ)  anode  serving  as  the  structural 
support,  a  YSZ  thin-film  electrolyte  layer  and  a  strontium-doped 
lanthanum  manganite  (LSM)-YSZ  cathode.  Although  there  are 
more  sophisticated  SOFC  designs  (e.g.  reactive  bi-layer  electrodes 
[32],  functionally  graded  electrodes  [33]),  we  choose  a  relatively 
simple  one  in  order  to  keep  the  degrees  of  freedom  in  design  as  low 
as  possible. 

The  cathode  electrode  consists  of  a  binary  mixture  of  electronic 
conducting  LSM  particles  and  ionic  conducting  YSZ  particles. 
Oxygen  (O2)  taken  from  atmospheric  air  is  transported  by  diffusion 
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dense  electrolyte  layer  (YSZ) 


from  the  cathode  gas  inlet  through  the  pore  phase  to  TPBs,  where 
the  oxygen  reduction  reaction  occurs  as 

^02  +  2e--01 2-.  (1) 

Electrons  are  provided  from  the  external  circuit,  and  the  oxygen 
ions  produced  are  transported  through  the  ionic  conducting  YSZ 
particle  networks  and  the  dense  YSZ  electrolyte  layer  towards  TPBs 
in  the  anode.  Besides  isolating  the  anode  and  cathode  for  electron 
conduction,  the  dense  electrolyte  layer  acts  as  a  gas  barrier.  In  the 
anode  cermet,  consisting  of  electron  conducting  Ni  particles  and 
ion  conducting  YSZ  particles,  the  hydrogen  oxidation  reaction  takes 
place  as 

H2  +  02  ^H20  +  2e-,  (2) 

where  hydrogen  (H2)  and  water  vapor  (H20)  are  transported  to 
and  from  the  active  TPBs  by  diffusion  through  pores,  and  the 
electrons  produced  are  conducted  through  the  Ni  particle  network 
to  the  current  collector  where  they  are  provided  to  the  external 
circuit. 


3.  The  extent  of  reaction  site  per  contact  pair  of  unlike  particles, 
given  where  the  electrochemical  reaction  takes  place  micro¬ 
scopically  at  the  conjunction  of  the  electron  conducting 
particle,  the  ion  conducting  particle  and  the  gas  phase. 

These  factors  will  be  discussed  in  detail  in  the  following. 

2.2. 1.  Average  coordination  number 

For  calculating  average  coordination  numbers  and  percolation 
probabilities  we  will  follow  closely  the  work  of  Chen  et  al.  [9]  and 
only  sketch  their  basic  ideas.  In  contrast  to  previous  micro-models 
for  estimating  the  coordination  number  of  binary  mixtures  [34,35], 
their  approach  also  satisfies  the  important  contact  number 
conservation  requirement  [5], 

In  a  binary  mixture,  the  surface-area  fractions  of  electron  con¬ 
ducting  particles,  Sei,  and  ion  conducting  particles,  Sj0,  are  related  to 
the  volume  fractions  c/>|(  and  the  particle  radii  r k  as 

5  j  hei4TCr2t  _  ^3) 

61  4tt  (nelr2  +nior?0)  <Pei/rei  +  q>'o/r''° 


2.2.  Effective  material  properties  in  anode  and  cathode 

We  consider  the  composite  anode  and  cathode  as  binary 
mixtures  of  spherical  particles.  Effective  microstructure  properties 
such  as  active  TPB  size,  effective  electronic  and  ionic  conductivity 
and  effective  pore  radii  for  the  gas  phase  can  be  determined  by 
using  micro-models  based  on  the  theory  of  particle  coordination 
for  randomly  packed  spheres  together  with  percolation  theory 
[5,9],  We  start  by  calculating  the  disposable  active  TPB  size  for 
electrochemical  reactions  of  a  specific  microstructure  that  depends 
on  three  major  factors: 


Si0  =  l-Sel,  (4) 

where  nei  and  ni0  are  the  number  of  electronic  and  ionic  particles 
per  unit  volume,  respectively.  Note  that  the  volume  fractions  are 
defined  relative  to  the  solid  material,  so  that  <pe \  +  <pi0  =  1 .  Then, 
the  coordination  numbers  among  particles  of  the  same  kind 
depend  on  the  surface-area  fractions  as 


=  Z0Sel 

(5) 

=  Z0Sio, 

(6) 

1.  The  coordination  number,  describing  the  average  number  of 
contacts  a  particular  particle  in  the  mixture  forms  with  its 
neighboring  particles.  Since  the  electrochemical  reaction  takes 
place  at  contacting  electronic  and  ionic  particles,  the  average 
coordination  number  among  particles  of  different  type  in 
a  binary  mixture  is  of  interest. 

2.  The  percolation  probabilities,  i.e.  the  likelihood  that  particles 
belong  to  cluster  networks  connecting  both  ends  of  the 
configurations,  thereby  forming  a  conductive  pathway.  This 
takes  into  consideration  that,  e.g.  in  an  anode,  TPBs  are  only 
active  for  the  electrochemical  reaction  if  the  electronic  particle 
involved  is  part  of  a  cluster  of  the  same  kind  of  particles  con¬ 
nected  to  the  current  collector,  and  the  ionic  particle  is  part  of 
an  ion  conducting  particle  network  forming  a  conductive 
pathway  to  the  dense  electrolyte  layer. 


where  Ze i_ei  describes  the  average  number  of  contacts  of  an  electron 
conducting  particle  with  neighboring  electron  conducting  particles, 
and  Z\0_ i0  the  average  number  of  contacts  of  an  ion  conducting 
particle  with  neighboring  ion  conducting  particles.  The  average 
coordination  number  of  all  particles  is  denoted  Zq  and  has  a  value  of 
6  for  a  binary  random  packing  of  spheres  [5,6,34],  Furthermore,  the 
coordination  numbers  among  particles  of  unlike  kind,  relevant  for 
calculating  the  effective  TPB  size,  can  be  estimated  as 


7  (l  +  (rel/rio)2y 


?  (l  +  (rio/rel)2y 


(7) 

(8) 
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Here,  Ze i_i0  is  the  average  number  of  contacts  an  electron  con¬ 
ducting  particle  has  with  ion  conducting  particles  in  its  vicinity, 
whereas  Zi0_e i  is  the  number  of  contacts  an  ion  conducting  particle 
has  with  electron  conducting  particles. 


contact  perimeter  line  of  overlapping  electronic  and  ionic  con¬ 
ducting  particles, 

“line  =  27ITc  =  27trincsin0inc.  (12) 


2.2.2.  Percolation  probability 

The  probability  that  an  electron  conducting  particle  in  the 
binary  mixture  belongs  to  a  network  of  similar  particles  connecting 
both  ends  of  the  configuration  can  be  given  as  [5,9] 

,,.(,-^1^)“)“  (9) 

where,  according  to  Kuo  and  Gupta  [36],  Zei_ei  =  1.764  represents 
the  percolation  threshold.  Below  this  threshold,  that  depends  both 
on  the  volume  fraction  and  particle  size  ratio  of  electronic  and  ionic 
particles,  no  electric  conductive  clusters  long  enough  to  connect 
both  ends  of  the  configuration  will  form.  Proceeding  in  the  same 
manner,  we  can  find  the  probability  for  a  continuous  percolating 
network  of  ion  conducting  particles  as 


(10) 


Since  various  recent  studies  [24-27]  indicate  that  the  electro¬ 
chemical  reaction  takes  place  on  the  surface  of  particles  within 
a  finite  width  from  the  contact  perimeter,  we  propose  a  new 
approach  for  calculating  the  extent  of  the  reactive  site  per  contact 
pair.  As  illustrated  in  Fig.  2b,  the  electrochemical  reaction  is  allowed 
to  take  place  on  spherical  caps  (SPC)  up  to  a  distance  ofxinc  from  the 
contact  perimeter  on  the  surface  of  the  inclusion  particle,  and  up  to 
a  distance  of  xmat  on  the  surface  of  the  matrix  particle: 

«spc  =  2TC^r?nccos0inc  -  r?nccos^^  +  0inc^ 

+  2tc  (r^costfmat  -  iiatcos  +  0mat)  )  (13) 

where  0mat  is  the  contact  angle  of  the  matrix  particle  and  given  as 


0mat  =  arcsin| 


csin0inc^ 


(14) 


2.2.3.  Extent  of  reactive  site  per  contact  pair 

The  main  difference  between  the  microstructural  percolation 
models  of  Costamagna  et  al.  [5]  and  Chen  et  al.  [9]  is  how  the  extent 
of  the  reactive  site  per  contact  pair  is  specified,  i.e.  where  the 
electrochemical  reaction  is  assumed  to  take  place  microscopically 
at  the  TPB.  Costamagna  et  al.  [5],  describing  the  TPB  as  an  area  per 
unit  volume  (m2/m3),  assume  that  the  reactive  site  is  the  contact 
area  between  the  electron  and  ion  conducting  particle.  For  the 
following  considerations,  it  is  important  to  clearly  define  the 
contact  angle  between  the  involved  particles.  Without  loss  of 
generality,  we  define  the  binary  mixture  of  electronic  and  ionic 
particles  to  consist  of  matrix  particles  (radius  rmat)  and  inclusion 
particles  of  smaller  size,  i.e.  rmc  <  rmat.  Then,  the  circle-like  contact 
area  between  two  particles  becomes 


2.2.4.  Effective  three-phase  boundary 

Finally,  the  effective  three-phase  boundary  can  be  calculated  as 
the  product  of  the  extent  of  reactive  sites  per  electronic  and  ionic 
particle  contact  pair  ax,  the  total  number  of  contacts  between 
electronic  and  ionic  particles  per  unit  volume,  Zei_j0  nei,  and  the 
likelihood  that  both  electronic  and  ionic  particles  belong  to 
percolated  clusters  Pei  Pj0  as 


Ax  =  axZel_ioneVPelPio 
4*  «x'Zio  el^io  Pe\P\o- 
where 


(15) 


nd  = 


(4/3)7cr|, 


and  ni0 


(i  -  e)<pfe 
(4/3)trr| 


(16) 


adrc  =  Ttr2  =  7r(rincsin0inc)2,  (11) 

where  rc  is  the  contact  radius  and  0inc  the  contact  angle  of  the 
smaller  particle  as  illustrated  in  Fig.  2a.  In  contrast,  Chen  et  al.  [9] 
formulate  the  TPB  in  their  percolation  model  as  a  length  per  unit 
volume  (m/m3)  as  earlier  introduced  by  Sunde  [14],  stating  that  the 
electrochemical  reaction  takes  place  on  the  (one-dimensional) 


represent  the  number  of  electronic  and  ionic  particles  per  unit 
volume  [5,7,9]  depending  on  the  porosity  e ,  i.e.  the  volume  fraction 
of  the  gas  phase.  As  can  be  seen,  the  contact  number  conservation 
requirement  Zei_i0nei  —  Zi0-eini0  is  fulfilled  by  Eq.  (15),  and  the 
extent  of  TPB  being  available  for  the  electrochemical  reactions  in  an 
electrode  microstructure  depends  on  the  volume  fractions,  the 
particle  radii  and  the  porosity.  The  subscript  x  in  the  definition  of  A 


matrix  particle  inclusion  particle 


b 


matrix  particle 


inclusion  particle 


Fig.  2.  Geometrical  illustration  of  the  different  three-pfmse  boundary  approaches. 
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denotes  which  approach  for  the  amount  of  reactive  sites  per 
contact  pair  is  chosen,  i.e.  whether  aCirc.  aime  or  aSpc  is  employed. 
Thus,  the  TPB  formulations  Acirc  and  Aspc  are  area  per  unit  volume 
(m2/m3)  parameters,  whereas  Aijne  is  a  length  per  unit  volume  (m/ 
m3)  quantity. 


2.2.5.  Effective  electronic  and  ionic  conductivity 

The  effective  electronic  conductivity  of  a  porous  electrode 
depends  on  its  porosity  e,  the  volume  fraction  cpel  of  electron  con¬ 
ducting  particles  and  the  percolation  probability  Pei,  i.e.  the  likeli¬ 
hood  that  electron  conducting  pathways  form.  It  can  be  estimated  as 

of  =  4r=  o°el((\-e)VelPelr  (17) 

Pe  1 

where  is  the  bulk  electronic  conductivity.  The  Bruggeman  factor  p 
takes  the  effect  of  tortuous  conduction  paths  into  account  and  is  set 
to  a  typical  value  of  g  =  1.5  in  this  study  [7,9,29],  Proceeding  in  the 
same  manner,  the  effective  ionic  conductivity  can  be  calculated  as 

of  =  ^  =  4(0  -  £)<PioPioY-  (18) 

2.2.6.  Effective  hydraulic  pore  radius 

In  order  to  include  the  effect  of  the  microstructure’s  morphology 
on  the  gas  diffusion  in  the  porous  electrode,  a  mean  hydraulic  pore 
radius  is  used.  According  to  Cannarozzo  et  al.  [  13  ]  it  can  be  written  as 


(19) 


and  depends  on  the  porosity,  the  volume  fractions  and  particle 
radii.  This  effective  pore  radius  will  be  employed  for  characterizing 
the  Knudsen  diffusion  when  we  model  the  gas  transport. 


2.3.  One-dimensional  cell  level  model 


In  order  to  systematically  investigate  the  impact  of  the 
morphological  parameters  for  the  anode  and  cathode  on  the 
performance  of  the  SOFC,  a  one-dimensional  cell  level  model  is 
used.  In  detail,  this  model  has  to  consider  the  following  physical 
phenomena: 

•  Mass  transport  of  the  gaseous  species  involved  through  the 
porous  anode  and  cathode:  diffusive  transport  of  reactants 
from  gas  channels  to  the  TPB  as  well  as  products  from  reaction 
sites  back  to  gas  channels. 

•  Electrochemical  reactions  occurring  at  the  TPB  reaction  sites. 

•  Conduction  of  electrons  from  reaction  sites  through  the  Ni- 
particle  network  to  the  current  collector  (anode)  and  from 
the  current  collector  through  the  LSM-particle  network  to  TPBs 
(cathode). 

•  Ionic  transport  of  oxide  ions  from  cathodic  TPBs  through  the 
respective  YSZ-partide  networks  and  the  dense  electrolyte 
layer  to  anodic  TPBs. 

•  Heat  generation  caused  by  ohmic,  activation  and  concentration 
losses  as  well  as  the  resulting  heat-transport. 

Basic  assumptions  of  the  model  are  uniform  total  gas  pressures 
in  the  porous  anode  and  cathode,  as  well  as  a  steady-state  opera¬ 
tion  of  the  fuel  cell,  i.e.  the  gas  flows  as  well  as  the  electronic  and 
ionic  currents  are  time-independent. 

2.3.1.  Governing  equations 

The  governing  equations  describing  the  migration  of  electrons 
and  ions  as  well  as  the  gas  transport  and  the  heat-transport  in 


anode  and  cathode  electrodes  take  the  form  of  Poisson  equations 
[6,13,37], 

V(-p5fV<?i°)  =-An!"  (20) 

V(-^ffV0el.n)  =  An/n  (21) 

VNa  =  ^  (22) 

v(-kfvr)  m  Qact.n  +  Qc onc.n  +  Qjoule.n,  (23) 

where  the  subscript  n  denotes  the  anode  or  cathode,  respectively. 
Further,  the  ionic  potential  the  electric  potential  <2>ei,  the  ohmic 
resistivities  peff,  the  local  Faradaic  current  densities  in,  the  rate  of 
mass  transport,  Na,  for  gas  species  a,  the  Faradaic  constant  F,  the 
absolute  temperature  T and  the  effective  thermal  conductivities  k„{f 
are  used.  The  right-hand  sides  represent  the  volumetric  source 
terms  and  describe  the  conversion  of  species  and  charge  associated 
with  the  electrochemical  reactions  as  well  as  the  generation  of  heat, 
Q,  due  to  loss  mechanisms  such  as  activation,  concentration  and 
ohmicjoule  losses.  For  the  dense  electrolyte  layer  separating  anode 
and  cathode,  all  source  terms  except  the  Joule  heating  vanish  and 
the  governing  transport  equations  reduce  to: 

(24) 

v(— Zcf^eVr)  =  Qjoulejyte-  (25) 

In  the  following,  the  modeling  of  the  electrochemical  processes, 
the  mass  transport  and  the  heat-transport  will  be  discussed  in 
detail. 


2.3.2.  Electrochemical  modeling 

The  local  Faradaic  current  densities  describe  the  charge  transfer 
rate  per  unit  TPB  between  the  electronic  and  ionic  current  corre¬ 
lated  with  the  electrochemical  redox  reactions  in  anode  and 
cathode.  They  can  be  written  in  the  form  of  Butler— Volmer  equa¬ 
tions  [6,10,38] 

ion  =  i°an  [exp (Vj^an)  -  exp ^  -  /Sf^an) ]  (26) 


ica  =  1'ca  [exp(j8pj|^ca)  “  exp(  -  »  (27) 

where  ft  and  0r  are  the  forward  and  reverse  charge  transfer  coef¬ 
ficients  in  anode  and  cathode,  R  is  the  universal  gas  constant  and 


are  the  reference  exchange  current  densities  [12,37,38]  having 
Arrhenius-type  relations  depending  on  the  local  partial  gas 
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pressures  of  hydrogen  (Ph2)  and  steam  (Ph2o)  in  the  anode  and 
oxygen  (po2 )  in  the  cathode.  The  reference  partial  gas  pressures  pref 
correspond  to  the  respective  gas  pressures  at  the  gas  inlets,  i.e.  at 
the  anode’s  or  cathode’s  surface.  Further,  EH2  and  E0l  are  the  acti¬ 
vation  energies  for  hydrogen  oxidation  and  oxygen  reduction, 
respectively.  The  parameters  yan  and  yca  are  used  as  fitting 
constants  in  the  model.  For  the  case  where  Eq.  (11)  or  (13)  are 
employed  for  calculating  A,  these  currents  in  Eqs.  (28)  and  (29)  are 
per  unit  TPB  area  (A  m~2),  whereas  they  are  currents  per  unit  TPB 
length  (A  m-1)  if  Eq.  (12)  is  used.  The  irreversible  activation  losses 
(also  commonly  termed  “polarization”  or  “overpotential”)  required 
to  overcome  the  reaction  barrier  energies  in  anode  and  cathode  are 
defined  as 


4an  =  ^-(*io  — *el,an) 

(30) 

4ca  =£aq-(<5el,ca  — ®io), 

(31) 

where 

„  1.06036 

~  /,  -tA  01561  4 

(S) 

1.76474 


7 


1.03587 

/l.52996kBT\ 


exp 


1  =  yftJp 


where  kB  is  the  Boltzmann  constant  and  T  and  ?  are  the  charac¬ 
teristic  Lennard-Jones  collision  diameter  and  potential  energies  of 
the  respective  gas  species  [40,42],  On  the  other  hand,  the  effective 
Knudsen  diffusion  coefficient  of  gas  species  a,  depending  on  the 
average  pore  radius  rpore,  can  be  estimated  as 


(38) 


Finally,  the  rate  of  mass  transport,  as  predicted  by  the  dusty-gas 
(32)  model  for  binary  gas  mixtures  [10,39],  becomes 


are  the  respective  equilibrium  potential  differences  [17]  at  the 
reaction  interfaces,  which  are  correlated  with  the  local  partial  gas 
pressures  and  the  Gibbs  free  energies  G  of  all  species  involved. 
Under  equilibrium  conditions  and  in  the  absence  of  current  flow 
the  SOFC’s  theoretical  open-circuit  potential  becomes  the  Nemst 
value  that  can  be  calculated  using  Eqs.  (32)  and  (33)  as 


ch2o  ~  Gh2  ~  0-5Gq2 


pff(p£f)° 

Ph2o 


2.3.3.  Modeling  of  mass  transport 

When  comparing  different  models  for  mass  transport  inside 
a  porous  SOFC  electrode,  Suwanwarangkul  et  al.  [39]  came  to  the 
conclusion  that  the  dusty-gas  model  is  suited  best  for  simulating 
the  gas  transport  of  a  binary  gas  mixture  through  a  porous 
medium.  It  is  assumed  that  the  spatial  gradient  of  the  total  gas 
pressure  in  the  porous  electrodes  is  insignificant,  so  that  the  total 
gas  pressures  of  the  anodic  H2-H2O  system,  p0an  =  Ph2  +Ph2o. 
and  the  cathodic  O2— N2  system,  p0an  =  Po2  +Pn2.  are  constant 
throughout  the  whole  depth  of  the  electrodes  and  only  diffusive 
gas  transport  phenomena  have  to  be  taken  into  account.  The 
dusty-gas  model  considers  two  diffusion  mechanisms:  molecular 
diffusion  is  dominant  for  large  pore  sizes  and  high  gas  pressures, 
whereas  Knudsen  diffusion  becomes  important  when  pore  sizes 
are  much  smaller  than  the  mean  free  path  of  the  gas  molecules. 
As  explicitly  shown  in  [40],  the  effective  molecular  diffusion 
coefficient  of  a  gas  mixture  containing  the  two  components  a  and 
is 


neff  3 

Daf>=m 


where  r  is  the  tortuosity  factor,  M  the  molecular  weight  and  Na 
Avogadro’s  number.  Furthermore,  the  collision  integral,  which  is 
based  on  a  Lennard-Jones  potential  [41  ]  becomes 


Nh  2  =  - 


1  -  (1  -  VMh2/MH2o)  -p- 

Po,an 


VPh2  (39) 


1 

i  -  (i  -  VMoJMNl)  ^ 

RT 

<N2 

Vpo2, 


(40) 


containing  both  molecular  and  Knudsen  diffusion  mechanisms. 
Note  that  due  to  our  previous  assumption  of  constant  total  gas 
pressures  in  the  porous  electrodes,  only  pH2  in  the  anode  and  p02  in 
the  cathode  have  to  be  considered  in  the  rate  of  mass  transport 
equations,  while  the  respective  remaining  partial  gas  pressures  can 
be  calculated  as  pH2o  =  Po.an  -  Ph2  and  pN2  =  p0,an  -  Po2  ■ 


2.3.4.  Thermal  modeling 

Most  SOFC  micro-models  found  in  literature  are  isothermal, 
assuming  a  constant  temperature  profile  throughout  the  whole 
configuration  [6,7,10,17,29],  while  only  a  few  explicitly  take  the 
generation  and  transport  of  heat  associated  with  various  loss 
mechanisms  into  account  [43],  In  this  study  we  consider  heat 
generation  due  to  activation  losses  (Qact),  concentration  losses 
(Qconc)  and  ohmic  losses,  causing  Joule  heating  Qj0uie-  The  volu¬ 
metric  heat  source  terms  of  the  porous  electrodes  can  be  written  as 
[43,44] 


Qact.n  =  PnAn'n  (41) 

Qconc, n  =  (C’ref  -  £nq)  An«n  (42) 

Qjoule,n=  (pfoffn)^V^io)2+(Pel!n)_1(V<Pelin)2,  (43) 

where  the  subscript  n  denotes  the  anode  or  cathode,  respectively. 
Eeq,|ef  represents  the  theoretical  open-circuit  potential  of  each  elec¬ 
trode,  adding  up  to  the  Nernst  potential  VNernst  =  E|3’ref  +  Eca’ref  as 
formulated  in  Eq.  (34).  For  the  dense  electrolyte  layer  only  Joule 
heating  due  to  the  ionic  current  has  to  be  taken  into  account  and  the 
heat  generation  terms  reduce  to 
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Qjoule.lyte  =  (pfojyte)  V«>io)2- 


(44) 


2.3.5.  Numerical  implementation 

In  order  to  numerically  solve  the  governing  equations  of  the 
one-dimensional  SOFC  model,  the  coupled  system  of  partial 
differential  Eqs.  (20)— (25)  was  implemented  in  the  commercial 
finite-element  software  COMSOL  Multiphysics  3.5a.  Representing 
the  degrees  of  freedom  in  the  finite-element  model,  <t> i0,  (J>e\,  Ph2. 
Po2  and  T  are  computed  in  the  procedure  that  is  used  to  integrate 
the  equations  and  thus  describe  the  solutions  of  the  boundary  value 
problems  considered.  The  respective  boundary  conditions  for  the 
anode  and  cathode  surface  are  listed  in  Table  1.  For  the  interior 
boundaries  continuity  of  flux  is  ensured. 

Throughout  this  work,  the  reference  temperature  is  Tref  =  1073  K 
on  the  exterior  boundaries.  Under  operation  conditions,  the  current 
density  of  the  cell  iop  is  specified  on  the  exterior  boundaries.  Then, 
since  the  electric  potential  $ei,an  is  set  to  zero  on  the  anode’s 
surface,  the  value  of  <Pei,ca  on  the  cathode’s  surface  corresponds  to 
the  steady-state  operating  cell  voltage  l/0p(z0p). 

In  order  to  parameterize  the  ID  SOFC  cell  model  we  use 
experimental  data  published  by  [28],  who  consider  a  laboratory 
anode-supported  SOFC  button  cell  operating  on  H2  fuel  diluted 
with  H20  vapor  and  02  taken  from  atmospheric  air.  All  model 
parameters  for  the  anode,  the  dense  electrolyte  layer  (DEL)  and  the 
cathode  were  chosen  to  correspond  to  the  reported  experimental 
values  as  much  as  possible  and  are  summarized  in  Table  2.  The 
mean  particle  diameters  and  volume  fractions,  which  were  not 
reported  in  the  original  work  of  Jiang  and  Virkar  [28],  are  estimated 
following  Jeon  et  al.  [7],  who  assume  a  monodisperse  mixture  of 
electron  and  ion  conducting  particles  of  equal  size  and  equal 
volume  fractions.  On  the  basis  of  these  model  parameters,  the 
effective  material  parameters  of  the  electrodes  are  calculated  as 
described  in  section  2.2.  For  the  effective  conductivities  of  Ni,  YSZ 
and  LSM  as  defined  in  Eqs.  (17)  and  (18),  the  temperature- 
dependent  bulk  conductivities  (in  Q-1  m  ')  are  taken  from  the 
literature  as  [29] 


ffvsz  =  ff°0 


°Ni  =  °el,an 


=  3.34  x  104  exp 


exp 

=  9.5  x  107 - 


(45) 


(46) 


°LSM  =  ^elxa 


(47) 


Table  1 

Exterior  boundary  conditions  us 
in  the  ID  cell  model. 


>n  the  anode  and  cathode  surface 


Boundary  conditions  or 


$eUn  =  0 
n-V<Pio  =  0 

=  PhI 


RTiop 


rp/:;f  *poz 

Boundary  conditions  on  cathode  surface 
"•«dlD  =  iop 
n-V<Pio  =  0 


(1  -  v/MH.i/MHao)(PH.,/Po,an 


-  P  -  V^MQiP^/POM) 


Table  2 

Parameters  used  in  the  ID  cell  model. 


Parameter 

Anode 

DEL 

Cathode 

Ref. 

L  [pm] 

1100 

20 

50 

[28,37] 

Po [Pa] 

101300 

101300 

[28] 

rei  [nm] 

500 

500 

[7,28] 

Ho  [nm] 

500 

500 

[7,28] 

■Pel 

0.5 

0.5 

[7,28] 

<Pio 

0.5 

0.5 

[7,28] 

0.4 

0.4 

[32,37] 

T 

5 

5 

[32,37] 

e*c  ["] 

15 

15 

[5] 

As  ft- 

0.5 

0.5 

[38] 

keB  [Wm  'r1] 

11 

2.7 

6 

[43] 

MHi  [mor’kg] 

2.016.10”3 

MHjo  [mor’kg] 

18.02.10-3 

M0j  [mor’kg] 

32.00.10-3 

MNj  [mor’kg] 

28.01. 10“3 

?h2  [K] 

59.7 

[10,42] 

?h2o  [K] 

809.1 

[10,42] 

fo,  [K] 

106.7 

[10,42] 

?n2  [K] 

71.4 

[10,42] 

Th2  [A] 

2.827 

[10,42] 

Th2o  [A] 

2.641 

[10,42] 

To2  [A] 

3.467 

[10,42] 

Tn2  [A] 

3.798 

[10,42] 

%  LI  mor’] 

120, 103 

[17] 

E02  U  mor’] 

130-103 

[17] 

Ch2  LI  mol  ’] 

-147.1  103 

[37] 

Gh2o  LI  mor’  ] 

-454.9 -103 

[37] 

G0i  [J  mol-’] 

-228.0- 103 

[37] 

C02  LI  mol  ’] 

-236.4- 103 

-236.4- 103 

[37] 

Xinc,  Xma,  [nm] 

29.2 

29.2 

Trirc.  Yspc  [A  m  ] 

1.13  105 

1.41 -104 

7 line  [Am  ’] 

7.31 -103 

9.12  -10-4 

Adrc.0  [nr1] 

8.38  -105 

8.38 -105 

Aline, 0  [m-2] 

1.30-1  O’2 

1.30-1 0’2 

ASpc.o  [rn  ’] 

8.38 -105 

8.38 -105 

By  utilizing  acjrc  from  Eq.  (11)  and  aiine  from  Eq.  (12),  we 
calculate  the  effective  extent  of  TPB  for  the  morphological 
parameters  listed  in  Table  2  and  defined  in  Eq.  (15),  both  as  an  area 
per  unit  volume  quantity  ACirc,o  (m2/m3)  and  a  length  per  unit 
volume  quantity  Aiine,o  (m/m3).  Then,  by  using  a  best-fitting 
approach,  we  adjust  the  parameters  of  our  model,  yCirc,an,  7drc,ca 
and  yiine.an.  Tiine.ca  in  Eqs.  (28)  and  (29)  to  the  experimental  data 
reported  by  Jiang  and  Virkar  [28].  Finally,  the  constants  Xjnc  and 
xmat  in  Eq.  (13),  defining  the  distance  from  the  contact  perimeter 
that  is  electrochemically  active  on  contacting  electronic  and  ionic 
particles  (see  section  2.2.3)  are  chosen  in  such  a  way  that  the 
reactive  area  per  contact  pair  aspc  is  in  agreement  with  acirc,  and 
consequently  Aspc  =  Acirc.  Fig.  3  provides  comparisons  between 
results  of  the  numerical  ID  SOFC  model  and  experimental  data 
reported  by  Jiang  and  Virkar  [28]  that  was  used  to  parameterize 
this  model. 

While  a  good  agreement  is  achieved  for  higher  current  densities, 
i.e.  in  the  operating  range  close  to  the  maximum  power  density, 
there  are  some  deviations  between  experiment  and  model  near  the 
open-circuit  voltage.  Similar  findings  have  been  reported  and  dis¬ 
cussed  by  Grew  et  al.  [37],  who  explain  the  deviations  by  possible 
fuel  leakage  in  the  experiment  and/or  pinholes  in  the  dense  elec¬ 
trolyte  layer.  Other  possible  contributions  to  the  deviation  might  lie 
in  the  nature  of  the  percolation  micro-model  that  assumes 
a  constant  coordination  number  of  6,  as  well  as  in  capturing  the 
electrochemical  reactions  by  use  only  of  the  Butler-Volmers 
equation.  Furthermore,  as  reported  by  Kenney  et  al.  [1],  the 
assumption  of  uniformly  sized,  monodisperse  electron  and  ion 
conducting  particles  in  both  electrodes  leads  to  higher  TPB  when 
compared  to  particle  mixtures  of  distributed  sizes  as  found  in  real 
electrodes. 
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3.  Results  and  discussion 

3.1,  Comparison  of  different  TPB  formulations 

In  a  first  step  we  compare  the  three  geometrical  approaches 
shown  in  section  2.2.3  for  calculating  the  extent  of  electrochemi- 
cally  reactive  sites  Acjrc,  Aiine  and  ASpC,  especially  regarding  their 
differing  scaling  behavior.  For  the  morphological  parameters 
rei  =  ri0  =  500  nm  and  <pe  1  =  <Pio  =  0.5  the  calculated  reference  values 
of  Ao  are  listed  at  the  bottom  of  Table  2.  The  values  found  for  Xjnc 
and  xmat  of  about  30  nm  are  in  agreement  with  findings  from 
studies  based  on  density  functional  theory  [24]  and  models 
coupling  elementary  chemistry,  electrochemical  oxidation  reaction 
kinetics  and  diffusion  mass  transport  phenomena  [26,27],  which 
indicate  that  depending  on  the  exact  reaction  mechanism  in  the 
electrochemistry,  the  reaction  takes  place  up  to  a  few  tens  of 
nanometers,  but  less  than  100  nm  [25]  from  the  contact  perimeter 
line  between  electron  and  ion  conducting  particles. 

For  a  mixture  of  electron  and  ion  conducting  particles  of  the 
same  size  (rei  —  r,0)  and  equal  volume  fraction  (<pei  —  <pi0  =  0.5), 
the  scaling  of  A  with  particle  size  is  shown  in  Fig.  4. 

In  order  to  ensure  comparability  between  the  different 
approaches,  the  TPB  is  normalized  by  the  respective  reference 
values  Ao  from  Table  2.  Significant  differences  in  the  scaling 
behavior  are  found:  ACirC  is  inversely  proportional  to  the  particle 
size,  Acirc  ^r-1,  while  Aijne  shows  an  inverse  quadratic  behavior, 
Aiine x  r_2-  In  contrast,  the  scaling  behavior  with  the  new  TPB 
approach,  Aspc,  coincides  with  Aune  for  relatively  big  particle  sizes, 
i.e.  when  the  particle  radii  rei,  n0  are  much  larger  than  the  elec- 
trochemically  active  domains  Xjnc,  xmat  on  the  particles’  surfaces, 
but  shows  an  increasing  deviation  for  smaller  particle  sizes,  where 

Aspect  b,b  >2. 

Of  special  interest  is  the  TPB’s  scaling  behavior  in  binary  particle 
mixtures  with  electron  and  ion  conducting  particles  of  different 
size.  We  define  the  particle  size  ratio 

R  =  ri0/rel  (48) 

as  the  ratio  of  the  radius  of  the  ion  conducting  particles  to  the 
radius  of  the  electron  conducting  particles.  We  then  select  for  study 
three  values  of  the  particle  size  ratio,  R,  namely  R  =  0.5, 1  &  4.  For 


(no  =  rei)  with  identical  volume  fractions  (?>ei  =  <Pia  =  0.5):  comparison  of  the  three- 
phase  boundary  measures  A^n.,  Aiine  and  Aspc. 

better  comparability  of  the  results  obtained  by  the  different 
approaches,  we  further  introduce  the  mean  particle  radius  of  the 
binary  mixture  as 

r  =  <Pi0r i0  +  <Peirei  (49) 

i.e.  the  volume  fraction  weighted  average  of  the  particle  radii. 
Fig.  5a  shows  the  normalized  extents  of  the  TPBs,  Adrc.  A|jne  and 
Aspc  from  Eq.  (15)  for  three  mixtures  of  different  particle  size  ratios. 
The  electron  and  ion  conducting  particle  radii  are  chosen  in  such 
a  way  that  the  mean  particle  radius  of  each  binary  mixture  is  the 
same  and  equal  to  500  nm  for  the  volume  fraction  of  the  electron 
conduction  particle  cp™3*  where  the  maximum  for  A  occurs.  As 
a  consequence,  the  particle  sizes  are  computed  as 

rei  =  t - \ - and  rio  =  RreI.  (50) 

(1  -  ?>Sax)R  +  Vet™ 
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where  r  =  500  nm,  and  <p™x  depends  on  the  particle  size  ratio  R. 
Fig.  5a  shows  the  results  as  a  function  of  <pe i,  the  volume  fraction  of 
electron  conducting  particles,  while  for  all  three  binary  particle 
mixtures  considered  the  particle  size  ratio  R  and  the  particle  radii 
n0  and  rei  remain  fixed.  Note  that  results  exist  only  for  limited 
ranges  of  the  volume  fraction,  rpei,  as  outside  those  ranges  either  the 
electronic  or  the  ionic  particles  fail  to  percolate.  The  case  R  =  1 
represents  the  reference  configuration  as  described  in  section  2.3.5, 
where  the  maximum  for  A  occurs  at  <p™x  =  <p^ax  =  0.5.  In  this 
case,  all  three  approaches  coincide,  thereby  calibrating  Xinc  and  xmat 
in  Eq.  (13).  For  binary  mixtures  with  different  particle  sizes,  i.e. 
R  =£  1,  the  maximum  in  A  for  all  three  approaches  shifts  away  from 
rpe i  —  0.5:  For  R  >  1,  i.e.  n0  is  larger  than  rei,  the  volume  fraction  of 
ion  conducting  particles  has  to  be  increased  and  that  for  electron 
conducting  particles  decreased  in  order  to  assure  percolation,  as 
previously  reported  [5,9,30],  The  resulting  maximum  for  A  for  R  =  4 
occurs  at  =  0.228  with  <p™ax  =  0.772,  and  the  particle  sizes  are 
re i  =  151  nm  and  ri0  =  604  nm.  When  R  <  1,  the  volume  fraction  of 
electronic  particles  must  be  increased  and  that  for  ionic  particles 
decreased  to  assure  percolation.  As  a  result,  the  maximum  for  A  for 
R  =  0.5  occurs  at  r/)])|lax  =  0.646  with  <p™x  =  0.354,  and  the  particle 
sizes  are  re i  =  608  nm  and  n0  =  304  nm.  A  significant  difference 
becomes  apparent  from  Fig.  5a  when  the  extent  of  the  TPB  is 
compared  for  three  mixtures  of  the  same  mean  particle  radius: 
While  the  maximum  achievable  values  of  Aijne  and  Aspc  show 
increases  for  R  =£  1  when  compared  to  the  reference  configuration 
with  R  —  1,  the  maximum  of  ACirc  for  R  ^  1  drops  below  its  refer¬ 
ence  value  of  R  =  1.  This  finding  is  of  importance,  since  it  indicates 
that  the  approach  of  using  Aiine  or  Aspc  will  predict  a  behavior 
contrary  to  that  found  when  ACjrc  is  employed  in  a  process  for 
determining  the  best  possible  size  ratio  between  electron  and  ion 
conducting  particles  in  an  anode  or  cathode. 

For  the  three  particle  size  ratios  considered,  Fig.  5b  illustrates 
the  relative  deviation  between  the  two  TPB  extents,  Aspc  and  Aiine, 
calculated  as 

A  _  (Aspc/Aspco)  —  (-A-line M-line.o)  cc-,i 

SPC''ine  '  (A|ine/A|ine,o)  (51) 

as  a  function  of  the  mean  particle  radius  r  of  the  binary  mixture.  For 
the  reference  configuration  (R  =  1,  r  =  500  nm)  Aspc  and  Aiine 


concur  as  discussed  before.  The  relative  deviation,  Aspc,iine,  caused 
by  taking  spherical  caps  on  the  particles  as  electrochemically  active 
surfaces  instead  of  confining  the  electrochemical  reactions  to 
a  one-dimensional  line,  amounts  to  6.5%  for  R  =  0.5  and  23%  for 
R  =  4  at  a  mean  particle  radius  of  F  =  500  nm.  Lowering  the  mean 
particle  radius  leads  to  a  (non-linear)  increase  of  the  relative 
deviations,  indicating  that  Aspc  grows  faster  than  Aiine. 

In  order  to  affirm  the  findings  that  ACirc  shows  considerably 
different  results  from  Aiine  and  Aspc  when  the  extent  of  the  TPB  is 
calculated  for  a  microstructure  of  particles  of  different  size,  Fig.  6 
shows  the  normalized  TPB  extent  as  a  function  of  the  volume 
fraction  <pet  and  the  particle  size  ratio  R  within  a  range  from  1  to  10. 

As  before,  the  mean  particle  radius  is  kept  constant  at 
r  =  500  nm,  so  that  according  to  Eqs.  (48)  and  (49)  the  radii  of  the 
electric  and  ionic  conducting  particles  become  dependent  on  R  and 
<Pe\  as 


(52) 


Due  to  the  fact  that  contact  number  conservation  is  fulfilled 
when  the  electron  and  ion  conducting  particles  are  interchanged 
(see  section  2.2.4),  A  is  unchanged  under  the  transformation  R — >-(1/ 
R)  and  ipe\— >(1  -  <pe |),  so  that  it  is  sufficient  to  consider  R  >  1.  Use  of 
Aspc  (Fig.  6a)  and  Aiine  (Fig.  6b)  for  the  TPB  extent  shows  similar 
behavior:  the  TPB  extent  increases  when  the  particle  size  ratio  R  is 
raised  and  the  mean  particle  radius  is  kept  constant.  However,  the 
results  diverge  increasingly  for  high  values  of  R:  Since  the  mean 
particle  size  is  fixed,  the  size  of  the  electric  conducting  particle 
decreases  for  an  increasing  R  and  vice  versa.  In  such  circumstances 
the  value  of  Aspc  takes  into  account  the  fact  that  an  increasing 
fraction  of  the  smaller  particle’s  surface  becomes  electrochemically 
active,  a  feature  that  is  neglected  in  the  case  of  Aijne,  causing  the 
deviation. 

In  contrast,  Acirc  (Fig.  6c)  exhibits  a  completely  different 
dependence  on  R  by  showing  a  maximum  for  a  mixture  of  equally 
sized  particles  with  the  same  volume  fraction.  Changing  the 
morphological  parameters  R  and  <pe\  in  this  case  always  results  in 
a  decrease  of  the  electrochemically  active  surface  when  compared 
to  the  reference  configuration.  A  saddle-point  value  of  Acirc  =  0.59 
Acirc,o  is  reached  at  R  -  3.25  for  cpel  just  greater  than  0.2.  Increasing 
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R  further  leads  to  a  slight  increase  of  ACirC,  but  even  for  unreal¬ 
istically  high  values  of  R  the  reference  value  of  ACirc,o  is  not 
achieved. 

As  described  in  section  2.2.3  the  calculation  of  ACirc  is  based  on 
the  assumption  that  the  electrochemically  active  surface  is 
proportional  to  the  direct  contact  area  between  electron  and  ion 
conducting  particles  that  is  not  accessible  to  the  gas  phase. 
Furthermore,  experiments  [21-23]  indicate  that  the  scaling 
behavior  of  the  TPB  extent  for  a  monodisperse  binary  particle 
mixture  follows  A^r~b,  b  ~  2,  instead  of  h  —  1  as  found  for  Acjrc. 
Hence,  we  suggest  that  the  approach  of  using  ACirc  for  the  TPB 
extent  is  not  well  suited  for  investigating  binary  mixtures  of 
particles  of  different  size.  While  the  approaches  of  using  Aiine  and 
Aspc  agree  with  the  experimentally  observed  scaling  behavior,  only 
Aspc  allows  for  the  electrochemical  reaction  to  take  place  on  the 
particles’  surfaces,  whereas  Aijne  hypothetically  confines  the 
reaction  to  the  one-dimensional  contact  perimeter  line.  Though 
both  approaches  coincide  for  binary  mixtures  of  large,  equally 
sized  particles,  this  difference  becomes  increasingly  important  for 
binary  mixtures  containing  particles  of  small  size,  i.e.  when  the 
electrochemical  reaction  takes  place  on  a  considerable  portion  of 
the  smaller  particle’s  surface.  Then,  the  approach  of  using  Aspc 


becomes  physically  more  sound  than  use  of  Aijne  and  we  suppose 
that  Aspc  predicts  the  TPB  extent  more  accurately.  As  a  conse¬ 
quence,  the  electrochemically  active  three-phase  boundary  surface 
will  be  calculated  from  AAspc  in  the  following,  unless  stated 
otherwise. 

3.2.  Parametric  study  of  morphological  parameters  in  anode  and 
cathode 

The  microstructure’s  morphology  in  anode  and  cathode  has 
a  crucial  impact  on  the  fuel  cell’s  performance.  In  the  following,  the 
influence  of  the  electronic  and  ionic  particle  size,  their  size  ratio 
and  their  volume  fractions  on  the  maximum  achievable  power 
density  is  investigated  and  ways  are  discussed  for  increasing  the 
performance  of  the  fuel  cell.  For  this,  the  one-dimensional  cell  level 
model  as  described  in  section  2.3  is  employed  with  all  parameters 
as  listed  in  Table  2.  Typical  operating  conditions  are  chosen  as 
boundary  conditions,  such  as  an  operating  temperature  of 
T  =  1073  K,  a  total  gas  pressure  of  105Pa  in  both  electrodes  and 
fractional  partial  gas  pressures  at  the  inlet  of  p[ f2f  =  0.95, 
p\f0  —  0.05  on  the  anode  side  (humified  hydrogen)  and 
p{f  =  0.20,  p[|f  =  0.80  on  the  cathode  side  (atmospheric  oxygen 


Fig.  6.  Normalized  three-phase  boundary  measures  Aspc,  Au„e  and  Adrc  depending  on  the  particle  size  ratio  R  and  volume  fraction  i/jel  for  binary  mixtures  of  constant  mean  particle 
radius  r  =  500  nm. 
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from  air).  This  reference  configuration  shows  a  maximum  achiev¬ 
able  power  density  of  Pmax  =  16,490  W  m-2  at  a  cell  voltage  of 
Vop  =  0.49  V  and  an  operating  current  density  of  iop  =  33,660  A  m-2. 

3.2.1.  Influence  of  particle  size  on  performance 

In  order  to  investigate  the  influence  of  the  particle  size  in  anode 
and  cathode,  binary  mixtures  of  electron  and  ion  conducting 
particles  of  equal  size  (R  =  1)  with  equal  volume  fractions 
<pel  =  ipl0  =  0.5  are  considered.  Results  are  shown  as  a  function  of 
particle  size  in  Fig.  7a  for  the  anode  and  in  Fig.  7b  for  the  cathode. 
Note  that  anode  and  cathode  are  examined  separately:  when 
properties  on  the  anode  side  are  considered,  the  microstructure 
parameters  of  the  cathode  are  kept  as  those  of  the  reference 
configuration  and  vice  versa;  recall  that  the  reference  configuration 
is  a  binary  mixtures  of  electron  and  ion  conducting  particles  with 
equal  volume  fractions  <pe\  —  <pi0  —  0.5  and  equal  particle  size 
(R  =  1 )  of  ri0  =  rei  =  500  nm.  It  can  be  seen  in  Fig.  7  that  the  power 
density  is  maximized  at  a  finite  particle  size.  On  the  anode  side, 
when  A  =  Aspc,  the  maximum  power  density  occurs  for  particles  of 
radius  re i  =  ri0  =  435  nm  and  is  =  16, 530  Wm-2,  a  small 
increase  of  0.2%  over  the  reference  configuration  in  which 
re i  =  ri0  =  500  nm  and  Pmax  =  16,490  W  m-2.  Monodisperse  anodes 
with  particle  radii  between  approximately  300  nm  and  500  nm 
show  almost  the  same  performance.  When  the  particle  size  in 
a  monodisperse  cathode  is  varied,  a  maximum  power  density  of 
Pmax  —  21 , 41 0  W  m-2  is  achieved  for  reirj0  =150  nm  for  A  =  Aspc, 
representing  an  increase  of  30%  compared  to  the  reference 
configuration  in  which  re i  =  n0  =  500  nm  and  the  power  density  is 
Pmax  =  16,490  W  nr2.  In  contrast  to  the  anode,  the  maximum 
occurs  for  relatively  small  particles  and  is  fairly  sharp.  The 
maximum  achievable  power  densities  P^“  and  P^Jiax  for  mono- 
disperse  electrodes  of  equal  volume  fraction  y>el  =  cpio  =  0.5  will 
be  used  as  reference  values  for  further  parametric  studies  on  the 
particle  size  ratio  R  and  the  volume  fractions  <pei&<pio.  Note  that  all 
results  are  given  and  compared  for  all  three  models  ACjrc,  Aijne  and 
Aspc  for  the  extent  of  the  TPB  in  Fig.  7a  and  b.  As  expected,  the 
results  for  both  Aiine  and  Aspc  are  qualitatively  similar,  whereas  the 
results  for  Acirc  deviate  significantly. 

3.2.2.  Influence  of  volume  fractions 

For  electrodes  consisting  of  a  binary  mixture  of  monodisperse 
particles,  the  volume  fractions  of  electron  and  ion  conducting 
particles  have  an  effect  on  the  cell’s  performance.  In  Fig.  8  the 


maximum  achievable  power  density  as  a  function  of  the  particle  size 
is  compared  for  two  cases:  in  the  reference  case,  the  volume  frac¬ 
tions  of  electron  and  ion  conducting  particles  are  chosen  to  be 
identical  so  that  <pel  =  <pio  =  0.5.  In  the  other  case,  the  volume 
fractions  <pe\  and  tpi0  —  1  <pel  are  varied  within  the  percolation 
limits  for  a  given  particle  size  until  the  largest  possible  power 
density  is  obtained.  All  computations  in  this  case  are  carried  out  with 
A  =  Aspc.  Fig.  8  shows  the  volume  fraction,  <pe i,  yielding  the  highest 
power  density,  as  well  as  the  resulting  value  of  the  power  density. 

Considering  the  anode  (Fig.  8a),  we  find  that  the  volume  fraction 
yielding  the  highest  power  density  is  approximately  <pe\  =  0.5  for 
small  particles  having  r  =  200  nm,  decreasing  to  <pe\  =  0.44  for 
r  =  1000  nm.  In  contrast,  for  the  cathode  electrode  (Fig.  8b) 
a  constant  volume  fraction  of  y>el  =  0.41 5  yields  the  highest  power 
density  within  the  range  of  r  =  200  nm  to  r  =  1000  nm,  dropping  to 
smaller  volume  fractions  only  for  particles  smaller  than  this  range. 
Since  the  volume  fractions  in  monodisperse  particle  mixtures 
influence  the  available  TPB  as  well  as  the  ionic  and  electronic 
conductivities  in  the  electrode,  but  not  the  pore  radius,  the  volume 
fractions  giving  the  maximum  power  density  represent  a  trade-off 
between  activation  and  Ohmic  losses,  and  depend  on  design 
parameters  (e.g.  thickness  of  the  electrode)  as  well  as  material 
parameters  (e.g.  bulk  conductivities  of  the  electronic  and  ionic 
conducting  phases).  For  our  specific  model,  the  power  density  can 
be  increased  by  1%  (r  =  1000  nm)  in  the  anode  and  by  4% 
(r  =  1000  nm)  in  the  cathode  when  compared  to  electrodes  with 
identical  volume  fractions.  When  compared  to  the  values  of 
and  P^“  found  in  section  3.2.1  for  the  particle  sizes  that  yield  the 
highest  power  density,  further  performance  gains  of  approx.  0.3% 
(r  =  435  nm,  <pet  =  0.46)  for  the  anode  and  2%  (r  =  150  nm, 
<pe i  =  0.415)  for  the  cathode  can  be  achieved  by  varying  the 
volume  fractions  within  the  percolation  limits.  The  benefit  of 
varying  the  volume  fractions  away  from  <pel  =  q>i0  =  0.5  are  thus 
seen  to  be  modest,  even  negligible. 

3.2.3.  Influence  of  size  ratio  and  volume  fractions 

Up  to  this  point,  only  monodisperse  particle  mixtures  with 
electron  and  ion  conducting  particles  of  equal  size  have  been 
considered.  However,  the  particle  size  ratio,  R,  between  electron 
and  ion  conducting  particles,  as  defined  in  Eq.  (48),  influences  also 
the  effective  material  properties  and  thus  the  cell’s  performance. 

In  order  to  quantify  this  influence,  Fig.  9  illustrates  the  effect  of 
differently  sized  electronic  and  ionic  particles  on  the  power  density. 


Fig.  7.  Influence  of  the  particle  size  in  monodisperse  particle  i 


1  (tio  =  rei)  with  identical  volume  fractions  ipio  =  <pd  =  0.5  on  the  maximum  achievable  power  density. 
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For  various  fixed  values  of  the  radius  of  the  ion  conducting  particle, 
n0,  the  maximum  achievable  power  density  is  shown  as  a  function  of 
the  radius  of  the  electron  conducting  particle  re  i.  For  that  purpose,  for 
every  combination  of  re\  and  r-l0  the  volume  fraction,  cpel,  which  yields 
the  highest  power  density,  is  determined  in  a  series  of  calculations. 
For  better  comparability,  the  maximum  power  densities  are 
normalized  by  P^“  and  (see  section  3.2.1)  and  results  for 
monodisperse  particles  (R  =  1)  with  equal  volume  fractions 
(<pei  —  <pio  =  0.5)  are  plotted  in  addition.  Such  calculations  are 
carried  out  for  both  anodes  and  cathodes  with  electronic  particle 
radii  ranging  between  rei  =  25  nm  and  re\  =  1000  nm,  with  the 
restriction  that  n0  >  25  nm,  and  volume  fractions  ranging  between 
^perc.iow  <.  ^  <  (-/,Perc>uPj  where  <p^erc’ ow  and  (p^erc,up  represent  the 
lower  and  upper  percolation  limit  for  electron  conducting  particles. 
This  study  shows  that  for  the  anode  (Fig.  9a)  as  well  as  for  the  cathode 
(Fig.  9b)  particle  mixtures  of  differently  sized  electric  and  ionic 
particles,  i.e.  R  =£  l.lead  to  increased  power  densities  when  compared 
to  mixtures  of  monodisperse  particles  with  equal  volume  fractions. 

Furthermore,  a  simultaneous  parametric  study  of  the  particle 
size  ratio  R  and  the  volume  fraction  <pe\  is  performed:  for  fixed  radii 


of  the  electron  conducting  particle  rei,  the  particle  size  ratio  R  and 
volume  fraction  <pel  are  varied  until  for  each  combination  of  re i,  R 
and  <pei  the  maximum  achievable  power  density  is  determined 
using  A  =  Aspc.  This  is  exemplified  by  Fig.  10  for  a  cathode  with 
re i  =  75  nm:  in  this  case  the  highest  power  density  is  achieved  for 
a  size  ratio  of  R  =  4.7  and  an  electric  particle  volume  fraction  of  <pA 
=  0.16,  resulting  in  a  power  density  of  23,610  W  m-2.  The  results  of 
the  simultaneous  parametric  study  on  particle  size  ratio  R  and 
volume  fraction  <pe f,  showing  their  influence  on  the  maximum 
achievable  power  density,  are  illustrated  in  Fig.  11.  As  a  function  of 
the  size  of  the  electron  conducting  particle,  it  shows  the  maximum 
power  density  that  can  be  obtained  when  the  particle  size  ratio  R 
and  volume  fraction  <pe\  are  varied  simultaneously.  In  this  study,  we 
impose  restrictions  on  the  minimum  particle  size  (rei,  n0  >  50  nm) 
and  the  particle  size  ratio  (1/10  <  R  <  10).  As  before,  the  maximum 
power  densities  are  normalized  by  and  P[|;“  and  results  for 
monodisperse  particles  (R  =  1)  with  equal  volume  fractions 
(<pel  =  c Pi0  =  0.5)  are  plotted  as  comparison. 

When  considering  the  anode  (Fig.  9a  and  11a),  we  find  that  for 
electron  conducting  particle  sizes  smaller  than  rei  =  450  nm 


Fig.  9.  For  various  fixed  values  of  the  ionic  particle  size,  ri0,  the  maximum  achievable  power  densities  in  anode  (a)  and  cathode  (b)  are  plotted  as  a  function  of  the  electronic  particle 
size,  rei.  For  each  combination  of  rei  and  ri0  the  volume  fraction  <pei  is  varied  within  the  percolation  limits  until  the  highest  power  density  is  identified.  The  power  densities  are 
normalized  by  and  as  determined  in  Fig.  7. 
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a  function  of  the  particle  size  ratio  R  and  the  volume  fraction  <pe\.  A  maximum  is 
reached  for  R  =  4.7  and  <pel  =  0.16,  indicated  by  the  spot. 

a  particle  size  ratio  of  R  >  1  becomes  favorable,  leading  to  signifi¬ 
cant  performance  increases  when  compared  to  the  reference  case 
of  monodisperse  particles.  However,  when  compared  to  P^“, 
arising  from  R  =  1  and  q>el  -  <pio  =  0.5  with  rei  =  ri0  =  435  nm,  only 
a  very  small  performance  gain  of  about  0.7%  (at  rei  =  250  nm, 
R  =  3.2,  <pel  =  0.27)  is  found.  For  electronic  particle  sizes  larger  than 
rei  =  450  nm,  the  simulations  indicate  that  particle  mixtures  of 
small  particle  size  ratio  (R  <  0.2)  become  favorable  and  lead  to  an 
increased  performance.  For  example,  for  a  mixture  of  rei  =  600  nm, 
ri0  =  60  nm  and  <pei  =  0.93,  a  performance  gain  of  3%  compared  to 
pniax  is  found.  This  performance  gain  becomes  even  larger  when 
smaller  ionic  particle  sizes  or  smaller  particle  size  ratios  are 
permitted. 


When  considering  the  cathode  (Fig.  9b  and  lib),  we  find  that  for 
electron  conducting  particle  sizes  smaller  than  rei  =  200  nm, 
a  particle  size  ratio  of  R  >  1  leads  to  an  increase  of  the  maximum 
achievable  power  density.  For  example,  for  a  binary  particle 
mixture  of  re i  —  50  nm,  ri0  =  420  nm  (R  =  8.4)  and  <pe \  —  0.1,  a  14.4% 
gain  of  the  power  density  is  found  compared  to  Pjj;“ ,  arising  from 
R  =  1  and  <pe\  —  ipi0  —  0.5  with  rei  =  rj0  =  150  nm.  This  perfor¬ 
mance  gain  becomes  even  larger  when  smaller  electronic  particle 
sizes  or  larger  particle  size  ratios  are  permitted.  In  a  size  range  for 
the  electronic  particle  between  200  nm  <  rei  <  250  nm,  the 
maximum  of  the  power  density  is  found  for  monodisperse  particle 
mixtures  (R  =  1 )  where  only  the  volume  fractions  are  varied  away 
from  <pel  =  q>io  =  0.5  (as  described  in  section  3.2.2).  Then,  for 
electronic  particle  sizes  larger  than  rei  =  250  nm,  binary  mixtures 
with  particle  size  ratios  of  R  <  1  become  favorable  when  compared 
to  the  monodisperse  reference  configuration  with  equal  particle 
sizes  and  volume  fractions.  In  contrast  to  our  results  for  the  anode, 
we  find  it  not  possible  to  increase  the  maximum  achievable  power 
density  compared  to  P°;“  for  electronic  particle  sizes  of 
rei  >  200  nm  when  we  simultaneously  vary  the  particle  size  ratio 
and  the  volume  fractions  within  the  described  ranges,  even  for 
particle  mixtures  of  very  small  particle  ratio  and  ionic  particle  radii 
as  small  as  n0  =  25  nm.  In  summary,  when  compared  to  mono¬ 
disperse  electrodes  with  electron  and  ion  conducting  particles  of 
same  size,  our  model  predicts  an  increased  performance  of  elec¬ 
trodes  with  electron  and  ion  conducting  particles  of  different  size 
(i.e.  R  ¥=  1 )  over  a  wide  range  of  particle  sizes. 

3.2.4.  Discussion  of  the  parametric  study 

In  order  to  interpret  the  results  of  the  parametric  study,  and  for 
better  understanding,  the  dependence  of  the  effective  material 
parameters  (three-phase  boundary  extent  A,  electric  conductivity 
ionic  conductivity  of0ff,  effective  pore  radius  rp0re)  upon  the 
particle  size  ratio  R  and  the  volume  fraction  cpei  for  a  fixed  electric 
particle  size  rei  is  plotted  qualitatively  in  Fig.  12.  Arrows  indicate  the 
direction  in  which  the  respective  parameters  increase. 

When  the  size  of  the  electron  conducting  particle  rei  is  kept 
constant,  the  extent  of  the  TPB  increases  for  small  particle  size 
ratios,  i.e.  rj0  <  rei,  and  high  volume  fractions,  <pe i,  ensuring  perco¬ 
lation,  a  trend  that  occurs  due  to  an  increase  in  the  number  of  ion 


Fig.  11.  Maximum  achievable  power  density  in  anode  (a)  and  cathode  (b)  for  a  simultaneous  variation  of  the  particle  size  ratio  R  and  the  volume  fraction  <pd  as  a  function  of  the 
electronic  particle  radius  rei.  Restrictions  on  the  minimum  particle  size  (rei,rj„  >  50  nm)  and  the  particle  size  ratio  (1/10  <  R  <  10)  are  imposed.  The  maximum  power  densities  are 
normalized  by  Pj(;“  and  P(J;“  as  determined  in  Fig.  7. 
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Fig.  12.  Material  parameters  qualitatively  depending  on  particle  size  ratio  R  and  volume  fraction  tpcl  for  re i  kept  constant:  three-phase  boundary  A,  conductivities  trei,  and  pore 
radius  r^.  Arrows  indicate  in  which  direction  the  respective  properties  increase. 


conducting  particles  and  therefore  of  the  electrochemically  active 
contacts.  The  direction  in  which  A  increases  in  Fig.  12a  is  therefore 
down  and  to  the  right  as  indicated  by  the  arrow.  In  contrast,  the 
pore  radius  increases  for  large  size  ratio,  R,  and  small  particle 
volume  fraction,  cpeh  leading  to  a  smaller  resistance  in  the  gas  flow 
and  therefore  to  smaller  concentration  losses  when  the  particle  size 
ratio  is  increased.  This  trend  is  apparent  in  Fig.  12d.  This  means  that 
when  the  size  ratio  is  changed  in  order  to  decrease  activation  losses 
by  increase  of  A,  the  concentration  losses  are  inevitably  increased 
and  vice  versa.  The  ion  and  electron  conductivities  mainly  depend 
on  the  volume  fractions  of  the  respective  phases,  so  that  an  increase 
of  one  conductivity  directly  causes  a  decrease  of  the  other,  char¬ 
acterized  by  the  opposite  trends  in  Fig.  12b  and  d.  Thus,  a  best 
possible  performance  is  only  achieved  when  the  overall  effect  of 
activation  losses,  Ohmic  losses  and  concentration  losses  is  mini¬ 
mized,  trading  off  the  disparate  trends  in  Fig.  12. 

As  presented  in  section  3.2.1,  our  model  utilizing  the  new  three- 
phase  boundary  approach  characterized  by  Aspc  predicts  that  for 
electrodes  consisting  of  a  monodisperse  (R  =  1)  binary  particle 
mixture  with  identical  volume  fractions  (<pe!  =  0.5),  a  maximum  in 
power  density  and  thus  in  performance  is  reached  for  a  well- 
defined  finite  particle  size,  which  is  in  agreement  with  previous 
studies  [6,13],  A  further  improvement  of  the  electrode’s  perfor¬ 
mance  can  be  achieved  by  subsequently  varying  the  volume  frac¬ 
tions  of  the  electron  and  ion  conducting  phases  within  the 
percolation  limits,  as  discussed  in  section  3.2.2,  tending  to  shift  the 
volume  fraction  of  the  electron  conducting  phase  to  values  below 
q>e i  —  0.5,  i.e.  increasing  the  volume  fraction  of  the  ion  conducting 
phase,  and  therefore  increasing  the  effective  ionic  conductivity. 
This  can  be  explained  by  the  difference  between  the  ionic  and 
electronic  bulk  conductivities,  which  involves  several  orders  of 
magnitude  (see  Eqs.  (45)— (47)).  Thus,  the  overall  Ohmic  losses  can 
be  reduced  by  increasing  the  ionic  conductivity  without  much 
penalty  arising  from  the  accompanying  decrease  in  the  electronic 
conductivity.  However,  a  decreasing  <pel  also  reduces  A,  so  that 
a  trade-off  between  activation  and  Ohmic  losses  must  be  found. 
Chan  et  al.  [30]  reported  similar  findings  from  use  of  a  cathode 
micro-model.  They  conclude  that  a  minimum  of  the  overall  losses 
and  thus  a  maximum  in  performance  can  be  achieved  for  R  =  1  and 
cpi0  -  0.6.  In  this  work,  we  show  that  the  ideal  volume  fraction 
leading  to  maximum  performance  is  not  independent  of  the 
particle  size  in  the  electrodes:  when  evaluating  the  functional 
dependency  of  the  volume  fraction  <pei  on  the  maximum  achievable 
power  density  in  an  anode,  we  find  that  the  volume  fraction  <pel  that 
maximizes  the  power  density  increases  with  decreasing  particle 
radius,  whereas  the  respective  volume  fraction  <pA  that  maximizes 
the  power  density  in  a  cathode  decreases  with  decreasing  particle 
radius.  This  opposite  trend  is  a  consequence  of  the  different 


thickness  of  anode  and  cathode  as  well  as  the  different  expansion  of 
the  reaction  zones  in  anode  and  cathode. 

Several  previous  theoretical  studies  on  the  impact  of  the  particle 
size  ratio  R  on  an  electrode’s  performance  stated  that  a  particle  size 
ratio  of  R  =  1,  i.e.  electron  and  ion  conducting  particles  of  the  same 
size,  should  be  chosen  in  order  to  achieve  the  maximum  perfor¬ 
mance  of  the  electrode  in  terms  of  power  density  [6,17,30,31]. 
Contrary  to  these  studies,  our  findings  shown  in  section  3.2.3 
indicate  that  the  performance  of  anode  and  cathode  can  be 
further  enhanced  by  varying  the  particle  size  ratio  away  from  R  =  1, 
and  enhancements  in  the  power  density  for  R  =£  1  are  found. 

When  the  performance  of  electrodes  that  consist  of  differently 
sized  particle  types,  i.e.  R  =£  1,  is  investigated,  the  approach  for 
calculating  the  effective  three-phase  boundary  is  crucial,  as  dis¬ 
cussed  earlier  in  section  3.1.  For  those  studies  employing  a  TPB 
approach  in  which  Adrc  characterizes  its  extent  [5,6],  a  particle  size 
ratio  R  =£  1  always  results  in  a  decrease  of  A,  so  that  such  studies 
will  predict  an  ideal  particle  size  ratio  of  R  =  1,  regardless  of  the 
particle  sizes  and  volume  fractions.  Other  studies  on  the  influence 
of  the  particle  size  ratio  used  a  TPB  approach  in  which  Aiine  char¬ 
acterizes  its  extent  [17,30,31],  and  exhibit  the  experimentally 
observed  scaling  behavior  (cf.  section  3.1). 

For  example,  in  a  parameter  study  for  a  cathode,  Chen  et  al.  [17] 
considered  a  constant  volume  fraction  </>e,  when  investigating  the 
impact  of  the  particle  size  ratio  R,  hence  finding  minimal  losses  in 
the  cathode  for  R  =  1.  However,  as  shown  in  our  work,  a  variation  of 
R  always  has  to  be  accomplished  by  a  simultaneous  variation  of  the 
volume  fractions  ye\  and  cp\0  —  1  -  ipe\  in  order  to  enhance  the 
power  density.  While  Chan  et  al.  [30]  and  Chen  et  al.  [31  ]  took  the 
influence  of  rpe]  into  account,  they  only  considered  a  few  specific 
combinations  of  re\  and  R,  from  which  they  concluded  that  an  ideal 
value  is  R  =  1.  Furthermore,  in  all  studies  previously  mentioned 
only  the  electrode’s  performance  for  a  specific  current  density  was 
considered. 

In  contrast,  we  have  systematically  investigated  the  perfor¬ 
mance  of  anode  and  cathode,  accounting  for  a  large  range  of 
particle  size  ratios,  volume  fractions  <pe\  and  particle  radii.  Our 
model,  employing  the  new  three-phase  boundary  approach  in 
which  Aspc  characterizes  its  extent,  predicts  that  the  performance 
of  composite  cathodes  can  be  enhanced  when  electron  conducting 
particles  with  rei  <  200  nm  and  bigger  ion  conducting  particles  are 
used,  i.e.  a  particle  size  ratio  R  >  1,  instead  of  monodisperse 
particles.  In  anodes,  significant  enhancements  are  only  possible  for 
rather  small  particle  size  ratios  (R  <  0.2). 

In  Fig.  13a  and  b,  cell  voltage  and  power  density  characteristics 
of  anodes  and  cathodes  are  plotted  as  functions  of  the  current 
density,  with  morphological  parameters  that  yield  highest  perfor¬ 
mance  in  the  examined  ranges  of  particle  sizes  and  volume 
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Fig.  13.  Cell  voltage  and  power  density  depending  on  the  current  density  for  anode  and  cathode  microstructures  with  morphological  parameters  that  yield  the  highest  power 
densities  as  identified  in  the  parametric  study.  Solid  line:  Parameters  yielding  the  maximum  power  density  when  particle  size  and  volume  fraction  are  modified  simultaneously. 
Dashed  line:  Parameters  yielding  the  maximum  power  density  when  only  the  volume  fractions  are  varied,  but  not  the  particle  size  ratio  so  that  R  =  1. 


fractions  as  identified  in  the  parametric  study.  Two  cases  are 
considered:  anodes  and  cathodes  with  parameters  yielding  the 
maximum  power  density  when  particle  size  ratio,  particle  size  and 
volume  fraction  are  varied  simultaneously  (solid  lines),  and  for 
comparison,  anodes  and  cathodes  with  parameters  yielding  the 
maximum  power  density  when  only  the  volume  fractions  are 
varied  (cf.  section  3.2.2),  but  not  the  particle  size  ratio  so  that  R  =  1 
(dashed  lines).  A  restriction  to  the  minimum  particle  size, 
rei, n0  >  50  nm,  as  well  as  to  the  particle  size  ratio,  0.1  <  R  <  10,  has 
been  imposed. 

For  the  anode,  our  simultaneous  parametric  study  of  particle 
size  ratio,  particle  size  and  volume  fraction  in  the  described  ranges 
identified  a  particle  size  ratio  R  =  0.1  and  a  volume  fraction  of  the 
electron  conducting  particles,  cpel  =  0.93,  as  the  morphological 
parameters  that  yield  the  highest  achievable  power  density.  This 
small  particle  size  ratio  of  R  =  0.1  and  the  high  volume  fraction  of 
electron  conducting  particles,  <pe\  =  0.93,  lead  to  a  reduced  ionic 
conductivity  and  a  smaller  mean  pore  radius  when  compared  to  the 
best  configuration  with  equally  sized  particles,  having 
Oo  =  rei  =  435  nm  (R  =  1 )  and  a  volume  fraction  of  the  electric  phase 
of  cpe |  —  0.46.  In  Fig.  13a  this  becomes  apparent  from  the  steeper 
decrease  of  both  the  cell  voltage  and  the  power  density  for  high 
current  densities,  and  a  reduced  maximum  current  density,  indi¬ 
cating  higher  concentration  losses.  This  is  compensated  by 
a  strongly  increased  extent  of  the  TPB  and  an  increased  electronic 
conductivity  (cf.  Fig.  12),  leading  to  an  overall  increase  of  the 
maximum  achievable  power  density.  For  the  anode-supported 
SOFC  design  considered,  variation  of  the  particle  size  ratio  away 
from  R  =  1  yields  enhancements  of  the  power  density  only  for 
particle  size  ratios  above  a  threshold  of  R  <  0.2,  while  within 
0.2  <  R  <  1  it  does  not  lead  to  an  improved  performance.  Note  that 
the  small  particle  size  ratio  necessary  might  already  approach  the 
limits  of  reliability  of  percolation  theory  [36,45],  so  that  we  suggest 
that  varying  the  particle  size  ratio  away  from  R  =  1  in  order  to 
increase  the  power  density  is  only  of  limited  value. 

For  the  cathode,  our  simultaneous  parametric  study  of  particle 
size  ratio,  particle  size  and  volume  fraction  in  the  described  ranges 
identified  a  particle  size  ratio  R  =  8.4  and  a  volume  fraction  of  the 
electron  conducting  particles  cpe\  =  0.1  as  the  morphological 
parameters  that  yield  the  highest  achievable  power  density.  This 
high  particle  size  ratio  of  R  =  8.4  and  the  low  volume  fraction  of  the 


electron  conducting  particles,  <pe\  =  0.1,  lead  to  an  increased  ionic 
conductivity,  a  larger  mean  pore  radius  and  an  increased  extent  of 
TPB  when  compared  to  the  configuration  with  equally  sized 
particles  of  rei  =  200  nm  and  a  volume  fraction  of  the  electric  phase 
of  cpe i  =  0.42,  where  only  the  particle  sizes  and  volume  fractions 
were  varied  in  the  respective  parametric  study  for  identifying  the 
highest  achievable  power  density,  while  the  particle  size  ratio  was 
kept  at  R  =  1.  In  Fig.  13b  this  becomes  apparent  from  the  increased 
maximum  current  density  and  from  the  decreased  slope  of  the  cell 
voltage,  indicating  reduced  activation,  concentration  and  Ohmic 
losses  when  compared  to  a  mixture  of  monodisperse  particles.  In 
contrast  to  the  aforementioned  anode,  our  parametric  study  indi¬ 
cates  that  a  variation  of  the  particle  size  ratio  away  from  R  =  1 
yields  enhancements  of  the  power  density  monotonically  for  R  >  1 
and  is  only  limited  by  the  size  of  the  smaller  particle  re i  in  the 
binary  mixture.  Hence,  we  suggest  that  for  certain  design  param¬ 
eters  the  particle  size  ratio  R  can  be  used  to  improve  the  perfor¬ 
mance  of  composite  cathodes,  when  the  volume  fractions  are 
varied  simultaneously. 

4.  Conclusion 

In  this  work  a  new  approach  for  calculating  the  effective  three- 
phase  boundary  extent  in  porous  anodes  and  cathodes  of  SOFCs  is 
proposed.  The  new  approach  is  motivated  by  the  assumption  that 
the  electrochemical  reaction  takes  place  on  the  surfaces  of  con¬ 
tacting  electron  and  ion  conducting  particles  within  a  certain  width 
of  the  contact  perimeter.  This  is  in  contrast  to  prior  approaches, 
which  either  assume  that  the  reactive  site  is  the  contact  area 
between  electron  and  ion  conducting  particles  or  confine  the 
electrochemical  reaction  to  the  one-dimensional  contact  perimeter 
line.  Hence,  the  new  approach  is  physically  more  sound.  The 
approach  introduced  is  compared  to  previous  ones  found  in  the 
literature,  and  its  advantages  for  calculating  effective  properties  of 
binary  particle  electrodes  with  particles  of  different  size  are  shown. 

Furthermore,  the  new  approach  for  calculating  the  effective 
three-phase  boundary  extent  is  employed  in  a  parametric  study, 
investigating  the  influence  of  the  morphological  parameters  of 
porous  anode  and  cathode  microstructures  on  their  performance. 
For  electrodes  consisting  of  a  binary  mixture  of  spherical  electron 
and  ion  conducting  particles,  the  influence  of  the  particle  size,  the 
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particle  size  ratio  and  the  volume  fractions  on  the  cell's  perfor¬ 
mance  is  investigated  by  use  of  a  one-dimensional  cell  level  model 
accounting  for  the  processes  of  gas  transport,  electron  and  ion 
conduction  as  well  as  the  electrochemical  reactions.  Based  on  this 
study,  we  suggest  that  under  certain  conditions  the  performance  of 
cathodes  can  be  enhanced  by  choosing  a  microstructure  with 
electronic  conducting  particles  of  size  different  from  that  of  the 
ionic  conducting  particles  when  compared  to  the  best  possible 
configuration  of  monodisperse  particles. 
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Nomenclature 

A  area  [m2] 

D  diffusivity  [m2  s”1] 

Eh2  activation  energy  for  H2  oxidation  [J  mol”1] 

Eq2  activation  energy  for  O2  reduction  [J  mol”1] 

Eeq  equilibrium  electric  potential  difference  at  reaction 
interface  [V] 

F  Faraday’s  constant,  96,485  [C  mol”1] 

Ca  Gibbs  free  energy  of  species  a  [J  mol”1] 
i  electrochemical  reaction  current  density  [A  m”2] 

i°  reference  exchange  current  density  [A  m”2] 

iop  operating  current  density  [A  m”2] 
k  thermal  conductivity  [W  m”1  K”1] 

L  length  [m] 

n k  number  of  k-particles  per  unit  volume 

Na  molar  flux  of  species  a  [mol  m”2  s”1] 

Na  Avogadro’s  constant,  6.02210”23  [mol”1] 

Ma  molecular  weight  of  species  a  [kg  mol”1  ] 
pa  partial  pressure  of  gas  species  a  [Pa] 

Po  total  pressure  [Pa] 

Pk  percolation  probability  of  k-particle 

rPore  mean  hydraulic  pore  radius  [m] 

rk  radius  of  k-particle  [m] 

R  diameter  ratio  of  particles 

R  universal  gas  constant,  8314  [J  mol”1  K”1] 

Q  heat  source  [W  m”2] 

T  temperature  [K] 

Vop  operating  cell  voltage  [V] 

z  moles  of  electrons  per  mole  of  reactant 

Z0  overall  average  coordination  number  of  all  solid  particles 

Zk  average  coordination  number  of  k-particle 

Zk_i  average  number  of  contacts  between  k-particles  and  1- 
particles 

Greek 

a  amount  of  reactive  sites  per  contact  pair,  [m]  or[m2] 

/?  charge  transfer  constant 

y  reference  exchange  current  factor 

e  porosity 

r 1  activation  overpotential  [V] 

Fa  collision  diameter  of  species  a  [m] 

A  effective  three-phase  boundary  (TPB)  in  [m  2]  or  [m1] 

ACirc  extent  of  TPB  based  on  approach  in  Eq.  (11)  [m”1] 

Aiine  extent  of  TPB  based  on  approach  in  Eq.  (12)  [m-2] 


Aspc  extent  of  TPB  based  on  approach  in  Eq.  (2.2.3)  [m1] 
fi  Bruggeman  factor 

<pk  volume  fraction  k-type  particles,  relative  to  solid  material 
<I>  electric  potential  [V] 

p  resistivity  [Q  m] 

a  conductivity  [Q”1  m”1] 

t  tortuosity 

8  contact  angle  between  electrode-  and  electrolyte 

particles 

Slap  collision  integral  for  species  a  and  (3 

Subscripts 
el  electronic 

io  ionic 

an  anode 

ca  cathode 

lyte  dense  electrolyte 

joule  joule  heating 

cone  concentration 

act  activation 

K  Knudsen 

Superscripts 

eff  effective 

ref  reference 

TPB  three-phase  boundary 


Appendix  A.  Supplementary  material 

Supplementary  material  associated  with  this  article  can  be 
found,  in  the  online  version,  at  doi:10.1016/j.jpowsour.2012.05.014. 
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